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SIGNAL-TO- NOISE RATIO ANALYSIS AND EVALUATION OF THE 


HADAMARD IMAGING TECHNIQUE 

Daniel J. Jobson, Stephen J. Katzberg, 
and Robert B. Spiers, Jr. 

Langley Research Center 


SUMMARY 

The Hadamard imaging technique, using binary coded masks to scan image 
space mechanically, enables large numbers of image elements to be sampled 
simultaneously, in contrast to the sequential sampling of conventional imaging. 
This technique was originally used to improve the signal-to-noise ratio per- 
formance of grating spectrometers. An analogous performance advantage when 
applied to imaging has not been verified. This paper analyzes the -signal- 
to-noise ratio performance of the Hadamard imaging technique and presents 
an experimental evaluation of a laboratory Hadamard imager. A comparison 
between the performances of Hadamard and conventional imaging techniques shows 
that the Hadamard technique is superior only when the imaging objective lens 
is required to have an effective F number of about 2 or slower. 


INTRODUCTION 

Electro-optical imagers used in most applications have measured the radi- 
ance distribution of a scene either by mechanically scanning an image across 
the active area of a single photodetector or by electronically scanning a 
stationary image with a detector array. In both types of imagers, elements of 
the image are sampled sequentially. An imaging technique which uses binary 
coded masks to sample large numbers of image elements simultaneously has been 
suggested. (See ref. 1.) 

This imaging technique uses masks with patterns of opaque and transparent 
elements arranged according to a predetermined binary code to scan a station- 
ary image mechanically. The motion of the mask causes different combinations 
of image elements to be transmitted to a photodetector which responds with an 
electrical signal proportional to the sum of the transmitted elements. An 
image is recovered from this series of electrical signals when the solution of 
the set of equations which mathematically represents the measurement process is 
computed . 

One binary code particularly useful as a mask pattern is a cyclical modi- 
fied Hadamard code. This code is selected for reasons of compactness, ease of 
data reduction, and unbiased efficient sampling. Cyclical codes eliminate the 
need for a unique mask pattern for each successive sample of the electronic • 
signal. The use of Hadamard codes allows simple fast recovery of image data 
from the electronic signals because the solution of the set of equations 
requires a simple matrix inversion. In addition, the use of Hadamard codes 



provides unbiased efficient sampling since each image element appears in about 
one-half of the summations. (See ref. 2.) 

Although the use of Hadamard masks in grating spectrometers provides 
improved signal-to-noise ratio performance in comparison with conventional 
scanning, there is no experimental verification that a Hadamard imager can 
achieve analogous gains in performance when compared with conventional imaging 
methods. In fact, one paper in the literature has expressed doubt that sub- 
stantial improvements can be achieved. (See ref. 3.) 

A general analysis of the signal-to-noise ratio performance of the 
Hadamard imaging technique is given here. In addition, an analysis and the 
results of performance measurements of a. laboratory Hadamard imager are pre- 
sented. The results of the analyses and measurements are used to draw a com- 
parison of the Hadamard and conventional imaging techniques. 


SYMBOLS 

Aj^j Hadamard code elements arranged to represent imaging process 

A 0 objective lens aperture area, m 2 

a detector area, m 2 

C,C' empirical constants for "1/f” noise 

c proportionality constant, v/m-Hz^ 2 

D lens aperture diameter or length of one side of square aperture, m 

d image diameter, m 

d' image diameter using immersion lens, m 

dp depth of photoconductor, m 

E electric field, V/m 

F focal length, m 

F# effective F number (F# = Z/D) V 

f frequency, Hz 

G electronic gain 

G' system electronic gain exclusive of photoconductor and preamplifier 

Gp photoconductor gain 

H e height of mask code element, m 
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h 

I 

i 

i* 


i NT 


K 

k 


Planck constant 

spectral irradiance on photoconductor, W/m 2 -ym 
current , A 

photo-induced current, A 
total noise current, A 
equivalent input noise voltage, V 
Boltzmann constant 
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L 
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N j 

5 J 


NER 


n 

n' 
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n t 

PR 


P 

q 

R 1 > R 2> R B 
R(X) 


Si ,s, 


system constant , 

W/m 2 -sr 

focused object distance, m 

focused lens to image distance, m 

length between photoconductor electrodes, m 

number of scene elements 

radiance of jth scene element, W/m 2 -sr 

spectral weighted average radiance of jth scene element, W/m 
noise equivalent radiance, W/m^-sr 
index of refraction 

volume density of photogenerated charge carriers, m - ^ 

noise voltage in ith electronic signal, V 

total volume density of charge carriers, m“ 2 

performance ratio of Hadamard imager to conventional imager 

number of samples 

charge of electron, C 

various circuit resistances, ohms 
photoconductor responsivity , V/W 
ith, jth electronic signal, V 


Sj electronic signal due to jth image element, V 

T temperature , K 

tp lifetime of charge carriers in photoconductor, sec 

Vgg bias voltage, V 

V Q output voltage, V 

V NA preamplifier equivalent input noise voltage, V/Hz 1 ^ 

Vfjf noise voltage due to feedback resistor, V/Hz 1 ^ 

vnx total preamplifier equivalent output noise voltage, v/ Hz - *^ 

V N1> V N2» V NB noise voltages due to various circuit resistors, V/Hz^^ 

W electronic bandwidth, Hz 

W e width of mask code element, m 

w width of photoconductor, m 

x distance from Hadamard mask to field stop along optical axis, m 

(see fig. 2) 

y distance from field stop to field lens measured along optical axis, m 

(see fig. 2) 

a(x) quantum efficiency of photoconductor, electron-hole pairs/photon 

e error signal, V 

X wavelength, um 

P mobility of charge carriers, m^/V-sec 

v optical frequency, Hz 

photoconductor resistance, ohms 

t wavelength independent transmissivity of optics 

t(x) spectral transmissivity of optics 

photons per second arriving at photoconductor per unit wavelength, 
(nm-sec) -1 

ft image element solid angle, sr 
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Subscripts: 

c conventional imager 

F field lens 

H Hadamard imager 

i,j,k integers 

o objective lens 

rms root mean square 

An estimate is denoted by < >. A bar over a symbol indicates the aver- 
age of the parameter . 


THEORETICAL COMPARISON OF CONVENTIONAL AND HADAMARD IMAGING TECHNIQUES 

The signal-to-noise ratio performance of both conventional and Hadamard 
imaging has been investigated to provide a basis for comparison between the 
two methods. Noise equivalent radiance was chosen as the figure of merit for 
the comparison. The dominant noise was assumed to be a detector noise which 
is signal and frequency independent . Performance parameters which may differ 
in the two cases are subscripted as c for conventional and H for Hadamard 


Conventional Imager 

For a conventional imager (see fig. 1), the electronic signal after ampli- 
fication is related to the radiance of the jth scene element by 

Sj = A Q PG J* 2 N j ( X ) t c (X) R ( X ) dX + nj (D 

where A 0 is the aperture area of the objective lens; ft, the solid angle sub- 
tended by a detector with respect to the objective lens; G, electronic gain; 
Nj(X), the spectral radiance elements of the scene; t c (X) , the spectral trans- 
missivity of the conventional imager optics; R(X), the spectral responsivity of 
the detector; and n^, the noise voltage occurring in the jth measurement of S. 

A spectrally weighted average radiance is defined as 


r 


(X) T C (X) R ( X ) dX 


N, = 


r-« 


(X) R ( X ) dX 


( 2 ) 
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L 



An estimate of this average is obtained from equations (1) and (2): 


S i - n i 

N,\ = — 

J/ GA 0 ni c 

where 

I c = J* 2 T C (A) R(A) dA 


( 3 ) 


The noise equivalent radiance or root-mean-square (rms) error in estimating N j 
from p measurements of each is 



(4a) 


or 

(NER >j,o * <“ 

where n. rmg □ c(a for most detectors, c being some proportionality 

constant, a c the detector area, and W the electronic bandwidth. 


Hadamard Imager 

For Hadamard imaging, a rotating mask is placed in the image plane of an 
objective lens. (See fig. 2.) A field stop defines the image field, and a 
field lens condenses the radiant power transmitted by the mask on to a detec- 
tor. This imaging process is represented in simplified form as 


s i = 


E # m n j 

j = 1 


(5a) 


without the presence of noise, or in matrix form 


[S] = k s [A][N] 


(5b) 
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where Sj_ is the column matrix of electronic signals from the detector, k s 
is a system constant, Nj is the column matrix of scene radiance elements, 
and m is the number of elements in the Hadamard code (as well as the upper 
limit for i and j). A i j is an m x m binary matrix whose first row j 
is the Hadamard code, each ’successive row being the previous row cycled by one 
code element. After m samples of have been taken, a set of m equa- 

tions in m unknowns has been formed . 

To recover the scene radiance distribution, the matrix equation (5b) is 
solved for [N] 


[N] = k s -1 [A] -1 [S] (6) 


where [A]“^ is a matrix inverse. 

The noise equivalent radiance for the Hadamard imager is determined by 
expanding equations (5b) and (6) to include: (1) the parameters which comprise 

the system response to radiance, (2) the spectral characteristics of the scene 
radiance and electro-optical system, and (3) the effect of detector noise. 

The inclusion of these factors in the summation form of equation (1) yields 

m pX 2 

S ± = A o nG 23 A i,j J Nj(X) t H (X) R(X) dX + n i (7a) 

j=1 ’ X 1 

where £2 is now the solid angle subtended by each mask code element. It is 
often more convenient in practice to use an equivalent form of equation (7a) 
such as 


m 

S i = 53 A i > j s j + n i (7b) 

j = 1 

where 

r x 2 

Sj = A q £2G J Nj(X) t H (X) R(X) dX 
M 

(that is, the signal contribution of the jth image element). Returning to the 

C X 2 

matrix form of this equation and letting % = ) Tg(X) R(X) dX and Nj be 

X 1 

defined as before 
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[S] = A o nGI H [A][N] + [n] 
An estimate of the radiance is then 



urVcs] - cni) 
V gi h 


( 8 ) 


and the error in this estimate is 


_N - <N>] 


[£|[n] 

V GI H 


(9) 


<1 

The root-mean-square value of [A]~ [n] 


2m 1/2 n 


i , rms 


m + 1 


and gives 


is derived in reference 3 as 


(NER) 


i,H 


2 m 1 /2 nj^ rms 

m + 1 A q ^GI h 


( 10 ) 


Comparison 

Because of the difference in the two types of imaging systems, two param- 
eters in the NER expressions may not be the same for the two cases. These 
are optical system transmissivity and detector area. For this comparison, the 
Hadamard imager optical system is assumed to differ from that of the conven- 
tional imager only by the addition of a field lens and the deletion of a scan 
mirror. In addition, the detector areas are likely to be different because for 
Hadamard imaging the detector area does not have to match the image element 
area. In fact, the use of a small detector area is difficult to realize in 
actual hardware. After these considerations have been accounted for, the per- 
formance ratio PR of the Hadamard imager to the conventional imager from 
equations (4b) and (10) (assuming (NER)^ = NER, t c (A) = T c > anci T H^ = X H^ 
becomes 


(NER) C V a c )1/2(m + 

PR = « (11a) 

< ner >H t c (a H ) 1/2 (2m 1/2 ) 


Because m is usually a large number for most imaging systems, a very close 
approximation is 
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PR 


T H (a c ) 


1 / 2 ffl 1/2 


T c (a H^ 


1 / 2 . 


(11b) 


Since the transmissivities of the optics of the Hadamard and conventional 
imager systems should not differ greatly, the ratio of the detector areas is 
of primary concern. Equation (11b) shows that the performance advantage gained 
by simultaneous sampling of image elements (the va}^/2 factor) is realized 
only to the extent that the Hadamard imager detector area is made to approach 
that of a conventional imager. This condition requires demagnifying the total 
image field to as small an area as possible. First, use of a field lens is 
assumed to condense radiation on to the detector. Then the minimum detector 
area, determined from the relationship of object to image size for the field 
lens (see fig. 2), yields a minimum image diameter of 


d o z F 


^ • Tr 


( 12 ) 


Since 


.' /2 d ( 

<V ,/2 = -5- 


and 


U u > ,/2 = 


« ,/2 »h 


for circular detectors, 


PR « 


T H d c m 


1/2 


2t c D o 1 f/ L F 


(13a) 


But the image field stop is also the field-lens aperture and gives D ? = 
and since the field lens images the objective lens on to the detector, 
Lp=Z 0 +x+y results. (See fig. 2.) Thus, 


d c m 


1/2 


PR « 


t h l o + x + y / l F 


2x, 


D„ 


(13b) 
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or with x + y « l Q and effective F number defined as F# □ Z/D, the per- 
formance ratio becomes 


PR 


IL 

2 ^c f# F 


(13c) 


This expression illustrates that relative performance of Hadamard imaging 
is not m^ n / 2 

in practice. Rather, t.he performance advantage is a much 
smaller value (F# 0 /F#p) because the Hadamard imager cannot usually be made to 
have a detector area as small as that of the conventional imager. 

Although the Hadamard imager optics include an immersion lens, the conven- 
tional imager can also use an immersion lens with the same relative reduction 
in detector area. This effect would then cancel in the performance ratio. 

The minimum possible detector area and highest performance ratio are achieved 
by using the lowest effective F number possible for the field lens. 


EXPERIMENTAL EVALUATION OF A HADAMARD IMAGER 
System Description 

The Hadamard imager which was used as a test instrument is shown in fig- 
ure 3- This instrument operates in the thermal infrared spectral range from 
8 to 14 pm. Figure 2 presents a schematic drawing of the germanium optical 
subsystem; figure 4 presents a block diagram of the electronics subsystem. 

Table I summarizes the imager design parameters. 

The optical subsystem consists of objective lens, the Hadamard encoding 
mask, a field lens, and an immersion lens mounted directly on the liquid- 
nitrogen cooled mercury cadmium telluride (HOT) detector. An image of the scene 
is formed by the objective lens on the encoding mask. The field stop limits the 
area of the encoding mask which can transmit radiation from image elements to 
the remainder of the optical subsystem. Since the field lens images the objec- 
tive lens on to the detector active area, an out-of-f'ocus image of the encoding 
mask superimposed on the image field is formed at the detector. Thus, the 
effect of the field lens is to mix the radiation from the transmitted image 
elements on to the detector active area. The immersion lens is used to reduce 
the required active area of the detector to minimize detector noise. An ac 
motor rotates the mask at 720 revolutions per minute (rpm) . This rotation 
requires the detector signal to be sampled in 25-usec intervals. A frame of 
image data is acquired in 25.6 msec. The Hadamard mask contains timing slots 
to synchronize detector signal sampling with mask rotation. These timing slots 
are used to create timing pulses by means of small tungsten lights and detec- 
tors on opposite sides of the mask. 

The electronics subsystem performs signal sampling, integration, amplifi- 
cation, and digitization. (See fig. 4.) In order to realize detector noise- 
limited performance, the rest of the electronics subsystem should introduce no 
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appreciable noise* Therefore, the analysis of the electronics subsystem centers 
on noise sources. 


Optical Subsystem 

Optics .- A detailed optical arrangement of the Hadamard transform imager 
is shown in figure 2. The objective lens of focal length F Q forms an image 
of some object on a digitally encoded field mask. The image frame area is 
defined by a field-stop aperture placed as close as practical behind the mask 
at some distance x. This square aperture transmits image radiation which 
corresponds to a mask area containing 31 by 33 elements. 

Upon passing through the field-stop aperture, the radiation is incident 
on a field lens placed a distance y immediately after the aperture. The 
best imager signal-to-noise ratio is realized when the field lens concentrates 
all the encoded radiation on a detector of the smallest size possible. The 
minimum image size the field -lens can form is proportional to the ratio of the 
image distance l p and the object lens distance Lp. Therefore, the smallest 
ratio will be obtained by using the shortest focal length or lowest F number 
lens. The lowest effective F number possible with refractive optics is 0.5. 
However, when refractive optics have F numbers less than 1.0, lens aberra- 
tions begin to negate any further size reduction. For this imager, an effec- 
tive F/1.38 lens was selected. 

The image diameter d is reduced further to allow a smaller detector 
size by using hemispherical immersion optics. When the detector is placed in 
optical contact with the flat back surface of a hemispherical lens of index of 
refraction n, the image diameter d can be reduced in size to 



Therefore, the higher the index of refraction of the immersion lens, the 
smaller d H can be made. In this imager, the immersion lens is germanium 
which has an index of 4 at a 10-um wavelength; therefore, the image diameter 
is reduced by 4. In order for the detector to be filled with radiation, the 
square detector diagonal w/2 is made equal to d H . 

Mask . - The coded mask consists of two arc sector areas placed diametri- 
cally opposite each other on a rotatable disk as shown in figure 5. The sec- 
tors are concentric with the rotation axis of the disk; the sector areas are 
centered and normal to the optical axis of the objective lens. Therefore, the 
sector areas are irradiated by the image alternately as the disk rotates. The 
sector areas are made of thin opaque metallic plates which have small elemen- 
tal areas chemically etched away in the form of a binary code. The transparent 
and opaque elemental areas are effectively square in shape with a width W e 
and height H Q . These elemental areas are arranged in 31 concentric rows. 

Each row contains 1088 elements, 65 of which are redundant. The other sector 
area is arranged identically. The modified Hadamard code given in reference 3 
determines which elements are opaque and which are transparent. Two image 
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frames of information are presented for processing with each rotation of the 
mask disk. The elemental solid angular field of view (in steradians) of the 
imager is given by 


a 


W e H e 



(15) 


Electronics Subsystem 

This section reviews the basic performance of the radiation detector, pre- 
sents an analysis of the succeeding preamplification and video processing, and 
compares those predictions with actual system measurements. 

Detecto r The detector used in this system is a commercial mercury cad- 
mium telluride (HCT) detector optimized for operation from approximately 8 to 
14 pm. Resistances of these detectors are a function of their temperature, 
field of view, and the object radiance. In the configuration used in this 
study the HCT is enclosed in a Dewar which maintains a detector temperature 
of 77 K and restricts detectable energy to a 60° field of view. Detector size 
is 2.5 mm by 2.5 mm, and it is directly bonded to a germanium immersion lens. 
The HCT used in this application operates in the photoconductive mode and has 
approximately 40 ohms of element resistance when viewing a 295-K field. 

Photoconductors such as the HCT very nearly approach background photon 
noise-limited performance at their peak wavelength. The basic photocon- 
ductor noise mechanisms are a "1/f" noise arising from contacts, generation- 
recombination noise (G-R), arising from random creation and removal of carriers 
from the conduction process, and finally, Johnson noise arising from the random 
motion of the free carriers themselves. (See ref. 4.) To put this discussion 
in more quantitative terms, the situation illustrated in figure 6 is used as 
the basis for a model of the photon detector and noise processes. The detector 
has a bias current impressed upon it and is assumed to be receiving a spectral 
irradiance I x in W/nT^-ym. The length between the electrodes is a and the 
photoconductor has depth d^ and width w. By assuming that the photocon- 
ductor absorbs the photons impinging upon it with quantum efficiency a s a(x), 
there will be an increase in equilibrium carrier concentration (per unit wave- 
length) of 


n' 


I x Jtwt p 


hvJtwdp 


(16a) 


n 1 = 

AW d p 


(16b) 
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where is the total number of photons arriving per second. The increase 

in current in the measuring circuit is then (ref. 5) 


q*, t aEpwd 

i* = (17a) 

Jiwdp 

G p 

i' = q$ x a — (17b) 


where Gp is Eyt and is called the photoconductor gain. The total current 
generated is then p 



In actual operation the photoconductor is generally operated with constant 
bias current. Thus an increase in photocurrent is exactly balanced by a 
decrease in voltage across the photoconductor to hold the current of the total 
device constant. Detection is therefore by measurement of the photoconductor 
voltage changes. Proper operation requires that these changes be small in 
comparison with the dark (nonsignal) bias current and voltage. Such a mode of 
operation is assumed here. 

As mentioned earlier, noise in photoconductors can be divided, into three 
major categories and may be represented as three equivalent noise current gen- 
erators per Hz 1 / 2 (see ref. 6): 




L 2 - 
l NT " 


4KTwqpnt<ip 


4(i’ + i 


bias 


)‘ 


c ’ (1 ' * iblas* 


/j* 2 o« x (JxVl + 2 .f(t p ) 2 J 


wdpfcf 


(19) 


where the first term on the right is Johnson noise, the second is G-R noise, 
and the third is 1/f noise. For the HCT detector used here, the basic 
signal-to-noise ratio resulting from dividing equation (18) by equation (19) 
approaches the value that would be expected for a perfect photoconductor. 


The figure of merit which best describes absolute detector perfor- 


mance D* takes on its maximum value greater than 1 * 10 


10 cm-Hz 1 ^ 2 per watt 
at about 12 pm for the HCT detector whereas an ideal photoconductor has a 
theoretical value of 3 x 10^ G cm-Hz^ 2 per watt at that wavelength. 
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Preamplifier.- In order not to degrade system performance from that 
possible with the HCT detector, the signal must be boosted by an extremely 
low noise preamplifier. This is a difficult requirement since the noise per 
Hz 1 ^ 2 presented to the preamplifier is an extremely low value. In fact, the 
voltage noise associated with the detector is not far from that of a resistor 
generating Johnson noise alone. Evaluating the Johnson noise component of 
equation (19) gives per Hz 1 ' 2 a value of 1.0 x 10“9 V for the 40-ohm HCT photo- 
conductor used in this system. If it is assumed that the G-R and 1/f compo- 
nents of noise are near this value, then the preamplifier noise must be less 
than 1 to 2 nV per Hz 1//2 . 


The preamplifier used in this system is a commercial low-noise operational 
amplifier whose noise figure is specified at 2 dB for a 100-ohm source imped- 
ance. Unfortunately, the bias on the detector requires a balancing circuit 
introducing additional noise. The circuit in figure 7 which represents the 
detector-preamplifier and balancing circuit can be analyzed to yield the con- 
dition for balance, sensitivity to signal variations, and noise. 

For low-frequency conditions and typical operating conditions, the per- 
formance of the circuit can be approximated by 


V o = ' V BB T. 


+ V 


BB 1 


R 2 + R l \ 

R 2 R 1 / + R B 


( 20 ) 


The condition for balance may be found by setting V Q □ 0 and solving for R 1 : 


R 


R B R 2 


( 21 ) 


Using the values for the circuit in this system (Rp = 50 ohms, Rg = 250 ohms, 
and = 40.5 ohms) yields a value for R^ of 308.6 ohms which is close to 
the measured value of 306 ohms. 


Since the preamplifier was supplied with an unknown feedback resistor R^., 
measurements were made to determine its value. Resistor R^ was varied, 
and the change in output signal was used in conjunction with equation (21) 
to give Rf = 4880 ohms. This result agrees well with the nominal value 
(5000 ohms ± 10 percent) given by the manufacturer. Note that equation (20) 
gives the sensitivity to signal changes; that is, 


9 V 


o 


9 ?X 


/ R 1 R 2 
\ R 1 + R 2 


( 22 ) 


where it is assumed that « Rg, as is the case. 



The effect at the output of the various noise sources can best be arrived 
at by combining the noise sources by using superposition and root sum of 
squares per Hz^ 1 2 



where 


+ 



(4kTR 2 ) 


1/2 


4kTR, 


R <1 1*2 


(23) 


In equation (19) 


qjjn^wd 


-, the detector noise approaches the Johnson 


noise value over the electrical frequencies of interest, and — ~ 1 + — . 


R2 R2 

The following equation thus represents the system output noise voltage 
per Hz^ 2 : 


'NT 


2 l Rb?X 


+ R- 


1/2 


(24) 


R B^X , 

Values for (= 34.48 ohms) and R 2 (= 43.0 ohms) yield (with 


r B + 


1 /P 

V NA = 0) a noise voltage per Hz ' of 


NT 


1 .27 x 10 -9 — 


(25) 


By assuming that the preamplifier noise figure of 2 dB is correct and that 
R^, □ 4880 ohms, one can expect the predicted preamplifier equivalent input 
noise voltage per Hz 1>/2 of 1.6 x 10“^ volt. 
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Table II gives measured equivalent input noise voltage obtained for the 

system with the HCT detector on the one hand and a 41 -ohm dummy resistor on the 

other hand. Note the good agreement of both cases with the above predicted 
value except at the lower frequencies. The influence of 1/f noise, which has 
been neglected, can easily explain the increase in noise at lower frequencies 
and is assumed to be the source of difference between low-frequency noise 
values and prediction. Since 1/f noise is an empirically determined quantity 
and not predictable, it is assumed to be defined by these measurements. 

The effect of 1/f noise can be taken into account by integrating the 
function C/f over the frequencies encountered. The value assigned to C 
is the observed value of 1/f noise at a frequency at which the 1/f is 

dominant. At a frequency of 5 kHz, the equivalent input noise voltage is 

4.3 * 10-9 per Hz or 

K / n \2 

— = (4.3 x 10-9) 

5 x 103 

K = 9.2 x 10- 11 * 


The noise voltage can be determined by integrating the squared noise 
voltage over the assumed frequency range between 40 Hz and 40 kHz, and taking 
the square root (assuming for the moment that 1/f is the only noise source). 
The complete effect of the two dominant noise sources, the white and 1/f com- 
ponents, can be computed with a root-sum-of-squares calculation. The value 
for the equivalent input 1/f noise is 0.799 * 10-6 y. The value for the 
equivalent input white noise component in the 40-kHz bandwidth is 0.400 yV. 
After these noise sources have been assembled, a root-sum-of-squares calcula- 
tion gives an equivalent input noise voltage of 0.900 * 10“6 V. 

Including the approximate preamplifier gain of 115 and a postamplifier 
gain of G' gives a predigitized noise of 1.035 x G* x 10” 1 * V. 


Signal and Noise Measurements 

A calibrated blackbody source positioned 1 . 02 m from the imager objective 
lens was set for 800 K and allowed to stabilize for at least 3 hours. A cir- 
cular aperture 0.87 cm in diameter was placed over the blackbody source. In 
addition, a known high-pass interference filter covered the aperture of the 
source to limit the lower wavelength of the radiation to 8 ym. The blackbody 
source aperture roughly corresponds in diameter to two image elements. Moving 
the source laterally and vertically positioned the image near the center of 
the image field to minimize the effect of off-axis aberrations and vignetting. 
Smaller lateral and vertical positioning was required to locate the image so 
that it filled one image element uniformly with roughly equal contributions to 
adjacent image elements. 
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A predicted value for the signal Sj for the central image element 
is given by equation (7b), with the values A Q = 9.133 * 10"^ m 2 , 

Q = 1.24 x 10~5 S r, G = 1206, Nj equivalent to the source radiance of 
an 800-K extended source, and t^U) the 0 P tics transmissivity. Calculations 
made by using a table of values for the wavelength dependence of a blackbody 
source radiance, the optical transmissivity, and detector responsivity give a 
value for the integral on the right-hand side of equation (7b) of 1.64 x 10^. 
Utilizing the electronic gain gives a predicted signal level for the central 
pixel of 0.224 V. 

An average signal voltage was determined for the central image element by 
taking the arithmetic mean of 50 measurements. This value was 0.24633 V with 
a standard deviation of 5.2 x 10"^ V. 

Observation of system operation indicated noise in excess of that expected 
from the detector-preamplifier-digitization process. In order to determine 
properly the effect of the Hadamard technique on system noise and to compare 
the noise measurements with prediction, three noise measurements were made. 
First, actual detector noise measurements after amplification and digitization 
were taken with the detector removed from the imager. The detector then 
viewed a 300-K stationary background with emissivity equivalent to that of the 
mask. Second, noise measurements were taken with the detector replaced in the 
imager; thus, the rotating mask and the capped lens are viewed. The final 
noise measurement was taken after processing through the Hadamard inversion by 
taking an average of the rms noise of 10 image elements. (Fifty measurements 
of each element were made.) The results are summarized in table III. Excess 
noise was being generated when the detector viewed the rotating mask. This 
noise was approximately a factor of three greater than the electronic noise, 
including quantization noise. 

The reduction in system noise was found to be a factor of 14.5 when com- 
pared with the system noise before Hadamard processing. Compared with the 
detector-amplifier-digitization noise, the reduction was only a factor of 5. 
Thus, the Hadamard process reduced system noise by a factor close to the pre- 
dicted value of 16 (that is, * m^ 2 /2 ). However, for the particular instrument 
used in these tests, the system noise was somewhat higher than it should have 
been. Investigation -indicated that vibration in the rather large code wheel 
was the source of this excess noise. Proper mechanical design should vir- 
tually eliminate this problem; therefore, it is not expected to be a basic 
limitation . 


COMPARISON OF EXPERIMENTAL AND THEORETICAL RESULTS 

The foregoing analysis and measurement results suggest that four specific 
signal-to-noise ratio performance values be compared. In addition, the analy- 
sis provides a means for generalizing the performance comparison between the- 
Hadamard imager and the conventional imager. 

First, the four specific signal-to-noise ratio performance values given 
for comparison are all for the case of one image element filled by an 800-K 
extended blackbody source; therefore, all measurements and predictions can be 
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directly compared. The first of these four values is measured performance 
for the laboratory Hadamard imager limited by excess noise; the second value 
is the measured performance of a Hadamard imager limited by detector noise • 
assumed to be reduced by the multiplex factor of 16. A prediction of the best 
possible Hadamard imager performance is given by the third value . This value , 
calculated by (1) using equation (7b) and table I to determine the signal and 
(2) reducing the measured detector noise by assuming a field lens with an 
effective F number of unity, could be used to minimize the detector area. 

The fourth value is the best achievable performance with the conventional 
imaging technique. This fourth value is calculated by dividing the third 
value by the performance ratio (eq. (13c)), an optical transmissivity of 0.64 
for the conventional imager and 0.51 for' the Hadamard imager and F# Q /F#p = 1 
being assumed. A discussion earlier in this paper pointed out that the immer- 
sion lens is omitted from the comparison since it could be used with both 
imaging techniques with the same reduction in detector area. 

Comparison of the four values given in table IV reveals that the labora- 
tory Hadamard imager could be improved substantially by eliminating excess 
noise and by using the minimum possible effective F number field lens. Such 
improvements would lead to performance approaching that of the conventional 
imager. This comparison also implies that at best the Hadamard imager can 
perform no better than a conventional imager when a minimum F number objec- 
tive lens is used. 

The comparison of the two imaging techniques is generalized by using the 
performance ratio equation repeated here for convenience. 


T H F# o 

2 t c F# F 


An examination of this equation discloses that the Hadamard imaging technique 
surpasses conventional methods only when the field lens can be a much lower 
effective F number than the objective lens to more than offset the factor 
of 2 and the usually slightly lower optical transmissivity of the Hadamard 
imager. Since the practical minimum F number of the field lens is about 1, 
a performance advantage is achieved when the objective lens is required to be 
an F number greater than about 2. For practical systems, this performance 
advantage is not usually large and does not usually approach the full advan- 
tage of m^ /<2 /2. For an F/5 imaging objective, an F/1 field lens, and 
T H = T c» ^ or exam Pl e > the performance advantage would be 2.5. 

This performance comparison has considered only signal- to-noise ratio. 

In practical applications, the two techniques would also be compared for dif- 
ferences in dynamic range requirements, data reduction complexity, mechanical 
complexity, and cost in determining the relative merits of the two imaging 
techniques. 
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CONCLUDING REMARKS 


The Hadamard imaging technique was described, and an analysis of its 
signal- to-noise ratio performance was presented. An expression was derived 
for comparing potential Hadamard performance with conventional imaging tech- 
niques. The optical and electronics subsystems of a laboratory Hadamard imager 
operating in the thermal infrared spectral region were analyzed, and the signal- 
to-noise ratio performance of the imager was measured. 

The general comparison of Hadamard and conventional imaging was based on 
a performance ratio which predicts a performance advantage proportional to the 
F number of the imaging objective lens, rather than the higher performance 
advantage usually given as one-half the square root of the number of image 
elements. This significantly lower advantage results because the practical 
minimum detector size (and noise) is much larger than that required for the 
conventional imager. In addition, the potential advantage is offset further 
by a factor of 2 associated with the encoding mask and slightly lower optical 
transmissivity for the Hadamard imager. Consequently, a performance advantage 
can be realized only when an F number greater than about 2 is required for 
the imaging objective lens. 

Performance measurements of a laboratory Hadamard imager with a fast 
objective lens (about F/1) confirmed predictions for the Hadamard imager. 
Although this particular imager was limited by excess noise caused by a vibra- 
tion of the rotating mask, the pure detector noise was measured separately to 
determine potential system performance. In addition, the signal produced by 
a blackbody source was measured to determine absolute system responsivity , and 
the noise reduction property of Hadamard imaging was verified. 

Results of the general comparison, supported by the laboratory measure- 
ments, imply that Hadamard imaging will be useful primarily when design con- 
siderations dictate the use of a slow F number objective lens for imaging. 


Langley Research Center 

National Aeronautics and Space Administration 
Hampton, VA 23665 
December 20, 1976 
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TABLE I.- HADAMARD IMAGER PARAMETERS 


Objective lens: 

Aperture diameter, D Q , cm 3.41 

Focal length, F 0 , cm 3.81 

In-focus object distance, L 0 , cm 102.20 

Image distance, 1 0 > cm 3.96 

Assumed transmissivity 0.85 

Encoding mask: 

Elemental solid angle, 0, sr 1.24 x 10”5 

Height and width of element, H e and W e , cm ' 1 - 39 x 10“^ 

Field stop: 

Height, cm 0.53^ 

Width, cm 0.568 

Angular vertical field, deg 7.26 

Angular horizontal field, deg 7.72 

Field lens: 

Aperture diameter, Dp, cm ‘ . 1.25 

Focal length, Fp, cm . . 1.25 

In-focus object distance, Lp, cm 4.72 

Image distance, ip, cm 1.72 

Assumed transmissivity 0.85 

Immersion lens: 

Index of refraction at 10 iim, n 4.005 

Detector (HCT): 

Length of one side, w, cm 0.25 

Peak responsivity of immersed detector at 11 pm, V/W 26.9 

Other: 

Mask to field stop distance, x, cm 0.25 

Field stop to field lens distance, y, cm 0.513 

Image diameter after immersion, dg, cm 0.307 

Assumed loss due to circular image overfilling square detector ... 0.82 

Peak transmissivity at 9 pm of 8.0-pm cutoff filter 0.88 
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TABLE II.- COMPARISON OF PREDICTION AND EXPERIMENT FOR PREAMPLIFIER 
INPUT NOISE VOLTAGE WITH V BB = 10.0 VOLTS 


Frequency , 

Predicted noise 

Measured noise voltage, v/kz^ 2 , of 

kHz 

voltage, 

— 

— 


V/Hzl7 2 

41 -ohm resistor 

HCT detector 

5 

1.60 x 10-9 

4.3 x 10-9 

4.3 x 10“ 9 

10 

1.60 

3.3 

3.3 

15 

1.60 

2.5 

3.0 

20 

1 .60 

2.5 

2.5 

25 

1.60 

2.15 

2.5 

30 

1.60 

2.15 

2.17 

35 

1.60 

2.15 

2. 17 

40 

1 .60 

2.15 

2.0 

45 

1.60 

1.6 

2.0 

50 

1.60 

1.6 

1.9 
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TABLE III.- NOISE MEASUREMENT RESULTS 


Detector-preamplifier noise 

before matrix inversion 

System noise before matrix inversion 

Noise after matrix inversion . . . 


Root-mean- square 
noise voltage, 
volts 


2.24 x io -3 
6.75 x 10-3 
4.62 x 10' 4 


TABLE IV.- PERFORMANCE COMPARISON 


Signal-to-noise ratio 
(800-K source) 


Measured for system-noise-limited 

Hadamard imager 530 

Measured for detector-noise-limited 

Hadamard imager 1600 

Predicted for Hadamard imager with 

minimum detector area 2200 


Predicted for conventional imager with 
equivalent objective lens and image 
element size 


5300 



Objective lens 


Detector array 


Scan mirror 




— Objective lens 


H 


Single detector 


(a) Detector array for electronic scanning. (b) Scan mirror /single detector combination 

for mechanical scanning. 


Figure 1.- Conventional imager. 




Figure 2.- Hadamard imager optical subsystem. 
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Figure 4.- Block diagram of electronics subsystem. 
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Figure 5.- Diagram of rotating disk containing mask code and timing slots. 
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